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Abstract 

In  this  paper,  we  study  a  nonautonomous  size  and  class  age  structured  epidemic  model  with 
nonlinear  and  nonlocal  boundary  conditions.  We  establish  a  comparison  principle  and  construct 
convergent  monotone  sequences  to  prove  the  existence  of  solutions.  Uniqueness  of  solutions  is 
also  established. 


1  Introduction 

In  previous  research  literature,  infectious  disease  dynamics  are  often  modeled  via  systems  of  ordi¬ 
nary  differential  equations  in  which  population  structures  such  as  age,  size,  sex,  etc.,  are  neglected. 
However,  there  are  many  cases  in  which  incorporating  some  of  the  population  structures  into  the 
model  may  provide  additional  and  important  information,  and  may  be  helpful  in  the  understanding 
of  disease  dynamics.  For  example,  it  was  suggested  by  the  author  in  [16]  that  the  age-structured 
seasonal  interaction  rates  of  measles  that  are  due  to  school  attendance  may  cause  oscillations  that 
are  not  described  in  corresponding  models  without  age  structure. 

In  this  paper,  we  consider  a  three  compartment  model  with  a  susceptible  state,  exposed  state 
(latent  state),  and  infectious  state  for  the  spread  of  a  directly  transmitted  disease  in  a  size-structured 
population.  We  assume  that  the  disease  under  consideration  is  fatal  so  that  no  recovery  occurs.  We 
also  assume  that  the  disease  is  transmitted  only  horizontally,  and  that  every  newborn  is  susceptible 
to  infection.  Note  that  for  many  diseases  there  may  be  some  delay  between  the  initial  latent  state 
and  infectious  state.  Furthermore,  it  is  not  biologically  realistic  to  expect  all  the  individuals  to 
progress  into  the  infectious  state  at  a  fixed  time  period  after  the  initial  latent  infection.  Additionally, 
individuals  in  the  infectious  state  may  have  dramatically  varying  mortality  rates  depending  on  the 
times  that  they  progress  into  the  infectious  state.  Hence,  in  addition  to  size  we  also  need  to  record 
the  class  age  or  residency  time  of  the  exposed  and  infectious  states,  i.e.,  the  length  of  time  the 
individuals  have  spent  in  their  present  state.  Thus,  the  individuals  in  the  exposed  and  infectious 
states  are  structured  by  two  internal  variables:  size  and  class  age. 

To  our  knowledge,  the  paper  [17]  by  Sinko  and  Streifer  was  the  first  appearing  in  the  literature 
on  structured  population  models  with  multiple  internal  variables  in  which  the  populations  are 
categorized  by  chronological  age  and  size.  Although  their  simplified  version  with  size  only  as  the 
internal  variable  (the  classical  size-structured  population  model)  has  been  studied  extensively  in 
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past  decades,  the  original  version  has  not  been  as  thoroughly  investigated.  We  have  found  structured 
models  with  multiple  internal  variables  for  describing  cell  dynamics  ([5,  10,  11,  12,  14])  and  epidemic 
dynamics  ([13,  18]).  For  example,  the  authors  in  [12]  studied  an  autonomous  linear  chronological 
age  and  size  structured  cell  population  model  that  has  both  normal  and  quiescent  individuals;  the 
theory  of  positive  operator  semigroups  was  used  to  show  that,  under  general  assumptions  about 
individual  behavior,  the  age-size  distribution  of  the  population  converges  to  a  stable  distribution. 
A  proliferating  cell  population  structured  by  chronological  age  and  maturity  was  studied  in  [10] 
where  existence  and  uniqueness  of  the  solutions  was  established  via  the  method  of  characteristics. 
This  model  was  then  extended  to  a  nonlinear  one  in  [11]  with  a  nonlinear  term  corresponding  to  loss 
of  individuals  due  to  crowding;  existence  and  asymptotic  behavior  of  solutions  were  studied  using 
semigroup  operator  theory.  The  authors  in  [5]  investigated  a  linear  population  model  describing 
cellular  division  with  both  quiescent  and  proliferating  cell  populations  structured  by  chronological 
age  and  maturity  of  the  cell  in  which  a  proliferating  phase  duration  depends  on  the  maturity  of 
cells.  The  existence  and  uniqueness  of  solutions  were  established  using  semigroup  theory  by  first 
summing  over  the  chronological  age  to  reduce  the  model  to  a  system  of  size-structured  delay  partial 
differential  equations  with  the  time  delay  depending  on  maturity. 

Modern  analysis  of  structured  epidemic  models  with  multiple  internal  variables  apparently  be¬ 
gan  with  the  work  by  Hoppensteadt  in  [13]  and  Waltman  in  [18],  where  the  population  is  structured 
by  chronological  age  and  class  age.  In  [15]  a  general  model  of  a  population  structured  by  several 
internal  variables  is  formulated  as  a  nonlinear  integral  equation,  and  existence,  uniqueness  and 
positivity  of  solutions  is  established  using  the  method  of  characteristics,  the  contraction  mapping 
theorem  and  semigroup  theory. 

The  goal  of  this  paper  is  to  establish  the  existence  and  uniqueness  of  the  solutions  to  the  system 
(2.1)  described  below  by  defining  upper  and  lower  solutions,  establishing  a  comparison  principle 
and  constructing  monotone  approximations.  Using  this  approach,  we  replace  the  true  solution  in 
all  the  nonlinear  and  nonlocal  terms  with  some  previous  estimate  for  the  solution.  We  then  solve 
the  resulting  linear  model  and  obtain  a  new  estimate  for  the  solution.  This  iterative  procedure 
yields  the  solution  of  the  original  problem  by  passing  to  the  limit.  The  method  also  provides  an 
explicit  solution  representation  for  each  of  the  iterates.  Thus,  an  efficient  numerical  scheme  can  be 
developed.  The  key  step  involved  is  a  comparison  principle  between  consecutive  estimates  [3].  This 
method  has  been  used  in  both  linear  and  nonlinear  size-structured  population  models  for  a  number 
of  systems  (e.g.,  [1,  2,  4]).  However,  we  believe  ours  is  the  first  paper  which  uses  this  method  to 
study  a  structured  model  with  multiple  internal  variables. 

The  models  we  consider  here  include  as  a  special  case  the  general  models  in  [6] .  These  models 
were  developed  as  part  of  a  biomass  vaccine  production  system  wherein  shrimp  protein  pathways 
are  co-opted  for  rapid  production  via  recombinant  protein  transfection.  In  [6]  the  authors  used 
simulation  studies  to  investigate  the  efficacy  of  such  a  production  system  as  a  first  response  coun¬ 
termeasure  to  toxic  attacks.  We  provide  here  a  well-posedness  theoretical  foundation  for  the  systems 
developed  in  these  applications. 

The  outline  of  this  paper  is  organized  as  follows:  We  present  our  model  in  Section  2.  We  then 
establish  a  comparison  principle  and  argue  the  uniqueness  of  solutions  of  our  model  in  Section  3. 
Monotone  sequences  are  constructed  in  Section  4  to  establish  existence  of  solutions  for  our  model. 
We  present  concluding  remarks  in  Section  5. 
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2  Size  and  Class  Age  Structured  Epidemic  Model 

The  model  we  consider  here  is  given  by 

r^max  /*00 

St(x,t)  +  (gs(x,t)S(x,t.))x  +  ms(x,t)S(x,t)  =  -S(x,t)  /  7  A (y,6)I(y,0,t)d0dy 

gS(0,t)S(0,t)  =  F(S(;t),E(;;t),I(;;t)) 

S(x,  0)  =  S°(x ) 

Et(x,  9,  t )  +  (gE(x,  t)E(x,  9 ,  i))x  +  Eg(x,  9 ,  t)  +  mE(x,  t)E(x ,  0,  t)  =  -pE(9)E(x,  9,  t ) 

r^max /*OC) 

E(x,0,t)  =  S(x,t)  J  J  X(y,  9)I(y,  9,  t)d0dy  C2-1) 

£(0, 0,  t)  =  0,  £(x,  0, 0)  =  £°(x,  0) 

0,  t)  +  {g\x,  t)I(x,  9,  t))x  +  Ig(x,  9,  t)  +  mI(9)I(x,  9,  t)  =  0 

COO 

l(x,0,t)  =  /  pE  (9)E(x,9,t)d9 
I(0,6,t)  =  o'}  I(x,9,0)  =  I°(x,9), 

where  (x,9,t)  €  [0,xmax]  x  [0,oo)  x  [0,T],  and 

P3C max  ^max  COO  r^max  /*00 

=  /  (3s (x)S(x,t)dx  +  /  /  PE(x,  9)E(x,  9,  t)d9dx  +  /  /  (31  (x,9)I(x,9,t)d9dx. 

Jo  Jo  Jo  Jo  Jo 

Here  S(x,t)  denotes  the  density  of  individuals  in  the  susceptible  state  having  size  x  at  time  t, 

E(x,  9 ,  t)  denotes  the  density  of  individuals  having  size  x  at  time  t  that  have  spent  0  time  units  in 

the  latent  (exposed)  state,  and  I(x,9,t)  denotes  the  density  of  individuals  having  size  x  at  time  t 
that  have  spent  9  time  units  in  the  infectious  state.  The  functions  gs(x,t),  gE(x,t)  and  gJ(x,  t) 
are  the  growth  rates  of  the  individuals  having  size  x  at  time  t  in  the  susceptible  state,  latent  state, 
and  infectious  state,  respectively.  The  functions  ms(x,t)  and  mE(x,t)  are  the  mortality  rates  of 
the  individuals  having  size  x  at  time  t  in  the  susceptible  state,  and  latent  state,  respectively,  and 
m1  (9)  is  the  mortality  rate  of  individuals  having  spent  9  time  units  in  the  infectious  state.  The 
function  /3s (x)  represents  the  fecundity  rate  of  individuals  having  size  x  in  the  susceptible  state. 
The  functions  f3E(x ,  9)  and  (31  (x,  9)  denote  the  fecundity  rates  of  the  individuals  having  size  x  that 
have  spent  9  time  units  in  the  latent  state,  and  infectious  state,  respectively.  The  function  A(x,  9) 
denotes  the  infection  rate  of  individuals  having  size  x  that  have  spent  9  time  units  in  the  infectious 
state.  The  function  pE(9)  represents  the  rate  at  which  individuals  having  spent  9  time  units  in  the 
latent  state  become  infectious. 

For  convenience,  we  will  use  the  following  notation  throughout  the  paper: 

Ds  =  (0,  xmax)  x  (0,  T),  Dei  =  (0,  xmax)  x  (0,  oo)  x  (0,T), 

DEg  =  (0,  xmax)  x  (0,  oo),  Dei  =  (0,  xmax)  x  (0,  n)  x  (0,T). 

The  space  Clr(DEI )  is  defined  by: 

{(j)  G  C1(Dei)  |  there  exists  some  constant  9p  >  0  such  that  </>(x,  9,t)  =  0  for  9  >  9p}, 
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and  the  space  Cqt(Dei)  is  given  by 

{(f>  £  Cl(DEI) |  there  exists  some  constant  0  <  9$  <  n  such  that  4>(x,  9,  t)  =  0  for  6  >  9$}. 

The  following  standing  hypotheses  will  be  assumed  throughout  the  paper: 

(HI)  The  individual  growth  rates  satisfy  g  ,gE  and  g1  £  Cl(Ds).  We  further  assume  that 
gS{x,t.)  >  0,  gE(x,  t)  >  0  and  g\x,t)  >  0  for  (x,t)  £  [0.  xmax)  x  [0,T]  with  gs(xmax,t)  =  0, 
gE(xma,Xlt)  =  0  and  g1  (xmax,t)  =  0  for  t  £  [0,T], 

(H2)  The  mortality  rates  ms ,mE  £  L°°(DS)  and  m1  £  L°°(0,oo)  are  nonnegative  functions. 

(H3)  The  fecundity  rates  /3s  £  L°°(0,  xmax)  and  /3E ,  (31  £  L°°((0,  x'ma,x)  x  (0,  oo))  are  all  nonnegative 
functions. 

(H4)  The  latent  to  infectious  rate  pE  £  L°°(0,oo)  is  a  nonnegative  function. 

(H5)  The  infection  rate  A  £  L°°((0,xmax)  x  (0,oo))  is  a  nonnegative  function. 

We  now  give  the  precise  definition  of  the  solution  for  the  model  (2.1). 

Definition  2.1.  (S(x,t),  E(x,9,t),  I(x,9,t))  is  a  solution  of  (2.1)  if  all  of  the  following  conditions 

hold: 

(i)  5  £  L°°((0,  T);  L1(0,  xmax)  fl  T°°(0,  xmax))  and  E,I  £  L°°((0,  T)]L1(DEq)  n  L°°(DEq)). 

(ii)  S(x,  0)  =  S°(x)  a.e.  in  (0,xmax),  E(x,9, 0)  =  E°(x,9)  and  I(x,9, 0)  =  I°(x,9)  a.e.  in 
(0,Xmax)  X  (0,  OO). 

(iii)  For  every  t  £  (0,T),  every  4>s  £  C|1(D‘S),  c f>E  £  C^r(DEI)  and  cj)1  £  C^r(DEI),  we  have 


S(x,  t)(f>s(x,  t)dx 


S(x,  0)cf)S(x,  0)dx  +  /  0s(O,s)F(S(-,s),E(-,-,s),I(-,-,s))ds 


t  rx  n 


+ 


<0±  Jo 

ft  fX  max 


[4>s(x,  s )  +  gS(x,  s)<j)x(x,  s )  —  ms(x,  s)cj)S (x,  s)]5(x,  s)dxd-s 


S(x,  s)(j)S (x,  s ) 


0  Jo 


Uo  Jo 


\{y,9)I{y,9,s)d9dy 


dxds. 


(2.2) 


f^rnax  POO 


E( x,  9 ,  t) 4>e(x,  9,  t)d9dx 


io  Jo 


f^max  /*00 


E(x,  9 ,  0)4>e(x,  9 , 0 )d9dx 


io  Jo 

(*t  /*3?max  r°° 


+ 


+ 


lo  J 0  JO 

f*t  /*#max  pOO 


[4> f(x,  9 ,  s)  +  gE(x,  s)(j)E ( x ,  9 ,  s)  +  4>E(x,  9 ,  s)]£'(x,  0,  s)d9dxds  (2.3) 


[mE(x,  s)  +  /oe(0)]£'(x,  0,  s)(j)E(x ,  0,  s)d9dxds 


lo  J 0  J 0 

ft  fX  max 


S(x,s)cf>  (x,  0,  s) 


o  J  o 


lyo  jo 


X{y,9)I(y,9,s)d9dy 


dxds , 


4 


and 


^max  POO 


I(x,  9 ,  t) 9 ,  t)d,Qdx 


o 

2-max  POO 


I(x,  9 ,  0 )4>I(x,  9 ,  0)d9dx  + 


t  fX  max  /*OC 

I  ( m  n 


/O  ,  J 0 


0  JO 


+ 


(j)  (x,0, s)  /  p  (9)E(x,9,s)d9dxds 

Jo  (2.4) 


'0,  j 0 

rt 


[<j)s(x,  6,  s)  +  g  (x,  s)<f>x(x,  0,  s )  +  ^(x,  6 ,  s)]/(x,  s)c?0  dxds 

: 

mI(9)I(x ,  0,  s)(/)I(x,  9,  s)d9dxds. 


>o  Jo 


3  Comparison  Principle  and  Uniqueness 

In  this  section,  we  give  the  definition  of  upper  and  lower  solution  of  (2.1)  and  establish  a  comparison 
principle.  Then  we  show  that  (2.1)  has  at  most  one  solution. 

Definition  3.1.  A  pair  of  functions  (S',  E,  I)  and  (S,  E,  /)  are  called  an  upper  solution  and  a  lower 
solution  of  (2.1),  respectively,  if  all  of  the  following  conditions  hold: 

(i)  S,S  G  T°°((0,  T);  T1(0,  xmax)  fl  L°°(0,  mmax)),  L°°((0,  T);  L\D%)  n  L°°(Df/)),  and 

(ii)  S(x,0)  <  S°(x),  S(x,  0)  >  S°(x)  a.e.  in  (0,xmax),  E(x,9,  0)  <  E°(x,9),  E(x,9,  0)  >  E°(x,9), 
L(x,9, 0)  <  I°(x,9 )  and  I(x,  0,0)  >  I°(x,6)  a.e.  in  (0,xmax)  x  (0,  oo). 

(iii)  For  every  t  G  (0,T),  every  nonnegative  function  t/’5  £  C^-D5),  V;E  £  Cq^zU7)  and  i/j1  G 
Cg^D®7),  we  have 


'o 


S(x,t)'i^  ( x,t)dx 


r  2-max  I"  t  _  _  _ 

>  /  S(x,  0)^5(a;,  0)cte  +  /  ^(O,  s)F(S(-,  s),  E(-,-,  s),  /(•,•,  s))c?s 

4o  ^  J  o 

pt  PX  max 

+  /  /  [V’f  (a:,  s)  +  <?5(x,  s)^f  (x,  s)  —  ms(x,  s)if)S(x,  s)]5(m,  s)dxds 

Jo  Jo 


t  PX  max 


S(m,  s)ips(x,  s) 


o  J  o 


X(y,9)I(y,9,s)d9dy 


o  Jo 


dxds. 


(3.1) 


S_(x,  t)'ips(x,  t)dx 


pxmax  pt 

<  /  S(x,  0)ips(x,  0)dx  +  /  'ips(0,s)F(S(-,s),E(-,-,s),I(-,-,s))ds 
Jo  Jo 

nXmax 

[ijjg(x,  s )  +  gs(x,  s)ij>x(x,  s)  —  ms(x,  s)ij>s(x ,  s)]A(m,  s)dxds 


t  rx„ 


S_(x,  s)ij)  (x,s) 


o  Jo 


OO 


X(y,  9)I(y,  9,  s)d9dy 


Uo  Jo 


dxds. 


(3.2) 
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E(x,  9,  t)ipE(x,  9 ,  t)d9dx 


E(x,  9, 0 )ipE(x,  9, 0 )d9dx 


1 0  Jo 

ft  fX max  /*00 


'0^0  Jo 

ft  fX max  fOO 


I0J0  Jo 

ft  fX  max 


[ijjg  ( x ,  Q,s)  +  g  (x,  s)^x  (x,  9 ,  5)  +  (x,  0,  s)]i?(x,  s)d6dxds 

[mE(x,  s )  +  pe(0)].E(x,  0,  s)ipE(x,  9,  s)d9dxds 


+  /  /  S(x,s)ip  (x,  0,  s)  /  /  A (y,9)I(y,6,s)d9dy  dxds. 

Jo  Jo  [Jo  Jo 


0  -'ll. 


E(x,  9,  t)ipE(x,  9,  t)d9dx 


E( x,  9, 0 )ipE(x,  9, 0 )d9dx 


1 0  Jo 

ft  fX  maxfOO 


'0  JO  Jo 

ft  /“^max  /*00 


[Vlf  (x,  9,  s )  +  gE(x,  s)ipE{x ,  0,  s)  +  ipE(x,  9,  s)]E(x ,  0,  s)d9dxd-s 

) 

[mE(x,  s)  +  p'E(6*)];B(.T,  0,  s)ipE(x,  9,  s)d9dxds 


io  Jo  Jo 

ft  fXm  ax 


+  /  /  S_{x,s)ip  (x,  0,  s)  /  /  A (y,9)I_(y,9,  s)d9dy  dxds. 

Jo  Jo  [Jo  Jo 


I(x,  9,  t)^1  (x,  9,  t)d9dx 


'o  Jo 

f^max  fOO 


I(x,  9,  0)ipI(x,  9 , 0 )d9dx 


/  o  Jo 

/*C  /*# max  /*00 


'0/(x,  0,  s) J  pE(0)E(x ,  0,  s)dOdxds 


to  Jo  Jo 

ft  fX  max  /‘OO 


[V>f(x,  9,  s )  +  g1  (x,  s)ip^.(x,  0,  s)  +  ipj)(x,  9,  s)]/(x,  0,  -s)d9dxds 

: 

m1  (9)I(x,  9,  s)ipI(x,  9 ,  s)d9dxds. 


o  '0  jo 


/(x,  0,  Pyip1 (x,  9 ,  t)d9dx 


o  Jo 


/  /  L{x,9,0)ip  (x,  9, 0)d9dx  +  /  /  ip  (x,0,s)  p  (9)E(x,  9,  s)d9dxds 

Jo  Jo  Jo  Jo  Jo  (o  fO 

rt  /■Xmax/'OO  V°-UJ 

+  /  /  /  9,  s)  +  g7(x,  s)ipx(x,  9,  s)  +  ipg(x,  9,  s)]/(x,  9,  s)d9dxds 

JoM Jo  Jo 

ft  fX  max  /‘OO 


/*^Cmax/*00 

/  /  /  mI(9)f(x, 9,  s)ipI(x,  9,  s)d9dxds. 

Jo  Jo  Jo 


Based  on  the  above  definition,  we  can  establish  the  following  comparison  principle. 


Theorem  3.2.  Suppose  that  all  the  hypotheses  (H1)-(H5)  hold.  Let  ( S,E,I )  and  ( S_,E,L )  a 
nonnegative  upper  solution  and  a  nonnegative  lower  solution  of  (2.1),  respectively.  Then  S  >  S_ 
a.e.  in  Ds ,  E  >  E  and  I  >  I  a.e.  in  DEI . 
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Proof.  Let  ua  =  S_  —  S,  =  E  —  E  and  u1  =  I  —  L  We  have  us(x,0 )  <  0  a.e.  in  [0,  xmax], 
uE(x,9,  0)  <  0  a.e.  in  [0,xmax]  x  [0,oo)  and  ^/(x,  9,0)  <  0  a.e.  in  [0,xmax]  x  [0,oo).  Choose 
nonnegative  functions  G  C'1(D'S),  G  C^r(DEI)  and  if1  G  C^r(DEI).  By  (3.1),  (3.2)  and 
us(x,  0)  <  0  a.e.  in  [0,  xmax],  we  find  that  for  every  t  G  (0 ,T) 


ub(x,t)ip  ( x,t)dx 

rt 


<  [  if>S(0 ,s)F(us(-,s),uE(-,-,s),uI(-,-,s))ds 
Jo  . 

ft  fXmax 


+ 


+ 


'0±  JO 

ft  fXmax 


[i/)g(x,  s )  +  gS(x,  s)ifx(x,  s )  —  ms(x,  s)ips(x,  s)]us(x ,  s)dxds  (3.7) 


/ o  Jo 

ft  fXma  x 


us(x,  s)ifs(x,  s) 
5(x,  s),ifs(x ,  s) 


2-max  /*00 


/o  J 0 

,E( 


/0  JO 

^max  /*CXD 


X{y,9)I(y,9,s)d9dy 
X{y, 9)uI(y, 9, s)d9dy 


dxds 

dxds. 


1 0  Jo 

By  (3.3),  (3.4)  and  u^(x,  0,0)  <  0  a.e.  in  [0,xmax]  x  [0,oo),  we  obtain  that  for  every  t  G  (0 ,  T) 


/o  Jo, 

< 


u^(x,9,t)ijj  ( x,9,t)d9dx 

t  /*3?max  /*00 

[Y>®(x,  0,  s)  +  gE(x ,  s)'ipE(x ,  #,  s)  +  ipE(x ,  0,  s)]tt'E(x,  0,  s)d9dxds 


'o  jo  J o 

ft  fXma x  /*00 


[mE(x,  s)  +  /oE(0)]w£'(x,  6*,  s)'il>E(x ,  0,  s)d9dxds 


10  J 0  JO 

/“^max 


(3.8) 


+ 

+ 


J  o 

r*c  /'Smax 


us(x,  s)'ipE(x ,  0,  s) 
5(x,  s)ipE(x,  0,  s) 


3?  max  /*00 


0  J  0 


U0  J  0 

^max  /*00 


X{y,9)I(y,9,s)d9dy 
X{y, 9)uI(y, 9, s)d9dy 


\_Jo  Jo 


dxds 

dxds. 


By  (3.5),  (3.6)  and  uT(x,  9, 0)  <  0  a.e.  in  [0,xmax]  x  [0,  oo),  we  have  that  for  every  t  G  (0 ,  T) 


3?  max  /*00 


'0  </0. 

< 


u1  (x,9,t)if  ( x,9,t)d9dx 

t  /*^max  /“OO 

-i/;7(x,0 ,  s)  /  pE(9)uE(x,9,s)d9dxds 


1 0  L0 

ft  fX max  /‘OO 


+ 


[V’f  (x,  0,  s)  +  g1  (x,  s)ifx(x,  9,  s )  +  iPq(x,  9,  s)]«7(x,  0,  s)d9dxds 


(3.9) 


/o  Jo  JO 

ft  /*2)max  /*00 


m1  (9)uJ (x,  9,  s)ipI(x,  9 ,  s)d9dxds. 


io  Jo  Jo 

Let  ifs(x,t)  =  exp (rst)(ps (x,t),  where  G  C,1(D5')  and  the  constant  ts  are  chosen  so  that 
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ts  >  ||m5||0O  +  || Alloc  sup 


3-max  /*00 


se[o,T]  l  Jo  Jo 


I(y,9,s)d9dy  >.  Then  by  (3.7),  we  have 


us(x,  t)ps{x ,  t)dx 


<  [  <PS{0,s)F(us(-,s),uE(-,-,s),uI(-,-,s))ds 

Jo  . 

pt  pXrnax. 


+ 

+ 

+ 


0J0 
t  rx n 


[(/?f  (x,  s )  +  gS(x ,  s)(pE(x,  s)]ir(x,  s)dxds 


(3.10) 


/  o  J  o 

pt  pX  max 


us(x,  s)(pS(x,  s) 
S(x,  s)(ps(x,  s ) 


n00  _ 

A(y,  0)/(y,  9,s)d9dy 

. 


dxds 


o  Jo 


2- max  POO 


A(y,  9)uT(y,  9,  s)d9dy 


UO  JO 


dxds. 


Let  ipE(x,9,t)  =  exp (rEt)ipE (x,9,t),  where  ipE  G  Cqt(Dei )  and  the  constant  te  are  chosen  so 
that  te  >  Ilm^Hoo  +  ||pE||oo.  Then  by  (3.8),  we  have 


'^max  pOO 


uE(x,  9 ,  t)<pE(x,  9,  t)d9dx 


I 0  J  0, 

pt  pXrnayi  POO 

< 


[<pf(x,  9,  s)  +  gE(x,  s)ipE (x,  9 ,  s)  +  <pE(x,  9,  s)\uE(x ,  9,  s)d9dxds 


/  o  To  J  o 

pt  pX max  /‘OO 


+ 

+ 

+ 


[te  —  ( mE(x ,  s)  +  pE(0))]itE(x,  0,  s)(pE(x,  0,  s)d9dxds 


>0  Jo  JO 

pt  pX  max 


(3.11) 


/()  J o 

pt  /*^max 


us(x,  s)<pE(x,  0,  s) 
S(x,  s)<pE(x,  0,  s ) 


2-max  POO 


UO  Jo 

2-max  POO 


A  (y,9)I(y,9,s)d9dy 
A  {y,9)uI{y,9,s)d9dy 


dxds 

dxds. 


>0  Jo  IJo  Jo 

Let  -07(x,0,t)  =  exp(ryt)<^7(x,  9,  t),  where  ip1  G  Cqt(Dei)  and  the  constant  r 1  are  chosen  so  that 
t 1  >  ||m^||oo-  Then  by  (3.9),  we  have 


max  POO 


u1  (x,  9,  Pjtp1  (x,  9,  t)d9dx 


o  Jo, 


t  rx„ 


< 


yr(x,0 ,s)  /  p  (9)u^(x,9,s)d9dxds 

1  o  Jp  Jo  (o  12') 

rt  flma-xfOO 

+  /  /  /  [ipl (x,  9,  s)  +  g\x,  s)<y9^(x,  9,  s)  +  <pg(x,  9,  s)]it7(:r,  9,  -s)d9dxds 

JoM  Jo  Jo 

P  t  PX  max  POO 

+  /  [r1  —  mI(9)]uI(x,  9,  s)pI(x,  9,  s)d9dxds. 

Jo  Jo  Jo 

We  choose  a  nonnegative  function  <ps  G  C,1(H5)  that  satisfies  the  following  problem 

(pE(x,s)  +  gs(x,s)(pE(x,s)  =  0,  0<s<t,  0<x<  xmax, 

[ P  (^maxj  s)  =  0,  0  <  s<t, 

(ps(x,t)  =  xS{x),  0  <  x  <  xmax, 


(3.13) 


where  xS  £  Cq°((0,  xmax))  and  0  <  \S  <  1-  We  choose  a  nonnegative  function  pE  G  Clr(DEI) 
that  satisfies  the  following  problem 


0 ,  s)  +  cjE  (x,  s)pE (x,  9,  s )  +  pE(x,  8,  s )  =  0,  (x,  0,  s)  <E  (0,  xmax)  x  (0,  n)  x  (0,  f) 
^E(Xmax,  9,  s)  =  0,  0  <  S  <  t,  0  <  9  <  71, 

<pE(x,  n,s)  =  0,  0  <  s  <  t,  0  <  x  <  xmax, 

<pE(x,  9,  t )  =  xE(x ,  0),  0  <  x  <  xmax,  0  <  9  <  n, 


where  xE  £  ^“((0,  xmax)  x  (0,  n))  and  0  <  xE  <  1.  We  choose  a  nonnegative  function  p1  G 
Cqt(Dei )  that  satisfies  the  following  problem 


pl(x,9,s)  +g\x,  s)pIx(x,9,s)  +  pIe(x,9,s )  =  0 ,(x,9,s)  G  (0,xmax)  x  (0,n)  x  (0,t), 
/(xMT,  0,  s)  =  0,  0  <  s  <  t,  0  <  0  <  n, 

<yj7(x,  n,  s)  =  0,  0  <  s  <  t,  0  <  x  <  xmax, 

p\x,9,t)  =  x \x,9),  0  <  x  <  xmax,  0  <  9  <  n, 


(3.15) 


where  x1  £  CTq0 ((0,  xmax)  x  (0,  n))  and  0  <  <  1.  The  existence  of  such  ips(x,s),  <pE(x,  9,  s), 

and  pJ(x,8,s)  follows  from  the  fact  that  all  the  above  problems  are  linear  with  local  boundary 
conditions.  The  boundary  and  initial  conditions  of  the  above  problems  imply  that  0  <  (ps,  ipE ,  p1  < 
1.  Substituting  (ps  that  satisfies  (3.13)  into  (3.10),  we  find 


lo 


us (x,  t)xS{x)dx 


t  rx  n 


t  ^max  POO 


<  pS’S  /  /  us+(x,  s)dxds  +  pS'E  II  I  uE+(x,9,  s)d9dxds 

Jo  Jo  Jo  Jo  Jo 

pt  P 3?max  P  OO 

+pSJ  II  /  uI+(x,9,  s)d9dxds, 

Jo  Jo  Jo 


(3.16) 


where  the  functions  us+,  uE+  and  uI+  denote  the  positive  parts  of  us ,  uE  and  u1 ,  respectively,  and 

S, 


pS’S  =  Halloo  '  ^ 


and  = 


•^max  POO  _ 

+  Tb  —  inf _ mb(x,t)  —  inf _ A (x,0)  inf  <f  /  /  I(y,  9,  s)d9dy  >,  pS'E  = 

(x,t)eDs  (a :,0)er>f/  se[o,T]  { J0  J0  J 

SJ  -  "qI"oo  +  II A||oo  sup  \  /  S(y,s)dy\$  where  Ds  =  [0,xmax]  x  [0,T]  and 
se[o,T]  Uo  J 

DEq  =  [0,xmax]  x  [0,oo).  Substituting  pE  that  satisfies  (3.14)  into  (3.11),  we  obtain 


r^max  Pf^ 


'0  J  0 

,E,S 


uE(x,  9,  t)xE(x,  9)d9dx 


t  rx  n 


t  PXma x  POO 


where  pE)S  = 
It 


<  p^’0  I  I  us+(x,s)dxds  +  pE’E  /  /  /  uE+(x,9,s)d9dxds 

Jo  Jo  Jo  Jo  Jo 

ft  Wmax/'OO 

+pE'J  /  /  /  uI+(x,9,  s)d9dxds, 

Jo  Jo  Jo 

P9C max  POO 


(3.17) 


loo  sup  <  /  /  I(y,  9,  s)d8dy  >,  pE’  =  r  -  inf  mE(x,t)  -  inf  p  (9) 

se[o,T]  Uo  io  J  ( x,t)eDs  6»e[o,oo) 

{P9C max  'j 

/  5(y,  s)dy>.  Substituting  ip1  that  satisfies  (3.15)  into  (3.12),  we 

Jo  J 
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have 


I  o  Jo 


<  I1 


I,E 


uI+(x,  9,  s)d9dxds, 


(3.18) 


u\x ,  0,  9)d0dx 

pt  pXma x /*00  /‘^max/*00 

/  uE+(x,0,  s)d9dxds  +  [i1,1  /  /  / 

Jo  Jo  Jo  jo 

with  the  constants  fiI,E  =  Hp^Hoo  and  /x /,y  =  t1  —  inf  m1  (9). 

0e[o,oo) 

Since  (3.16)  holds  for  every  x  ,  we  can  choose  a  sequence  {x  ,fc}fcLi  converging  to  %  ,  where 


1  if  us(x,t)  >  0 

0  otherwise 


XS(x)  = 

Substituting  xS’k  into  (3.16)  and  using  the  dominated  convergence  theorem,  we  find 


us+(x ,  t)dx 


<  d 


s,s 


t  PXrx 


us+(x,  s)dxd-s  +  dS'E 


t  roo 


+d‘ 


S,I 


/o  Jo 

/"Imax  rOO 


0  Jo  Jo 


uE+(x,9,  s)d9dxds  (3.19) 


uI+(x,  9,  s)d9dxds. 


o  Jo  Jo 


Similarly  we  find  (3.17)  holds  for  every  xE ■  Hence,  we  can  choose  a  sequence  {xE,k}kLi  converging 
to  xE i  where 

1  if  uE{ x,  9,  t)  >  0 

0  otherwise 


XE(x,9)  = 


Substituting  such  xE’k  into  (3-17),  using  the  dominated  convergence  theorem,  and  using  the  fact 
that  the  constants  fiE'S ,  ^E,E  and  [iE)I  are  independent  of  n,  we  find 


C^max  /*00 


uE+(x,  9,  t)d9dx 


'o  Jo 


<  d 


E,S 


t  rxn 


us+(x,  s)dxds  +  dE'E 


t  rx max  roo 


0  Jo 


+d 


E,I 


t  rx max  roo 


uE+(x,9,s)d9dxds  (3.20) 


o  Jo  Jo 


uI+(x ,  9 ,  s)d9dxds. 


o  Jo  Jo 


Similarly,  we  can  find 


f^max  dOO 


10  JO 


<  d 


I,E 


uI+(x,  9,  t)d9dx 

''t  r^max  TOO  rt  /“iCmax /“OO 

/  uE+(x,9,  s)d9dxds  +  [i1'1  '  '  '  ',,/+ 


o  Jo  Jo 


(3.21) 


u1  +  (x,  9,  s)d9dxds. 


o  Jo  Jo 


roo 

By  the  fact  /  e~dd9  =  1  and  (3. 19)— (3.21) ,  we  have 

Jo 


noo 

[e~eus+{x1 1)  +  uE+(x ,  #,  f)  +  -u/+(x,  0,  t)]d0dx 

/*c  rX max  POO 

<p  /  [e~eus+(x,  s)  +  uE+(x,9,s)  +  uI+(x,9,  s)]d9dxds, 

Jo  Jo  Jo 
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where  n  =  max{^s,s  +  gE,S ,  p S,E  +  pE,E  +  ftI,E ,  nS}I  +  pEJ  +  ft1'1}.  By  Gronwall’s  inequality,  we 
have 

^max  POO 

/  /  [e~eus+(x,  t)  +  uE+(x,  9,  t)  +  uI+(x,  8,  t)]dOdx  =  0, 


10  Jo 

,.S  s'  n 


which  implies  that  tr  <  0  a.e.  in  Z)  ,  u  <  0  and  u1  <  0  a.e.  in  .D^  .  □ 

Remark  3.3.  From  the  proof  of  Theorem  3.2,  it  is  easily  seen  that  for  any  function  w  £ 
L°°((0,  T);  L1(0,  xmax)  fl  L°°(0,  xmax))  with  ws(x,0)  <  0  a.e.  in  [0,  xmax],  if  the  following  inequality 
holds  for  every  t  £  (0,  T)  and  every  nonnegative  function  T5  £  C1(I?'S) 


/  wb(x,  t)^fb(x,  t)dx  <  /  / 

/o  io  io 


['hf  (x,  s)  +  g(x,  (x,  s)]rcs'(x,  s)dxds 

f^max 

As(x,  s)ws(x ,  s)\ks'(x,  s)dxd-s, 


(3.22) 


with  ,4s  £  L°°(DS ),  then  ws(x,t )  <  0  a.e.  in  H5. 

Similarly,  for  any  function  wE  £  L°°((0,  T);  L1((0,  xmax)  x  (0,  oo))  D  -L°°((0,  xmax)  x  (0,  oo))) 
with  wE(x,6,  0)  <  0  a.e.  in  [0,xmax]  x  [0,oo),  if  the  following  inequality  holds  for  every  t  £  (0,T) 
and  every  nonnegative  function  £  Cqt(Dei) 

f*X max  /*00 

/  /  wE(x,9,t)'S>E(x,9,t)d9dx 


t  /“^max  /‘OO 


<  /  /  /  [\Hf  (x,  9,  s)  +  g(x,  0,  s)  +  'hf  (x,  9,  s)]wE(x,  0,  s)d9dxds  (3.23) 

io  io  io 

!*t  /*3?max/*00 

+  /  /  /  AE(x,9,  s)wE (x,9,  s)'S>E(x,9,  s)d0dxds, 

io  io  io 

with  £lE  £  L°°(Dei),  then  wE(x,8,t)  <  0  a.e.  in  H'5'1. 

This  remark  will  be  used  in  Section  4.  Next  we  will  show  that  model  (2.1)  has  at  most  one 
nonnegative  solution. 

Theorem  3.4.  Suppose  that  the  assumptions  (H1)-(H5)  hold.  Then  there  exists  at  most  one 
nonnegative  solution  for  the  system  of  equations  (2.1). 

Proof.  Suppose  that  there  exist  two  nonnegative  solutions  (Si,  E\,  I\)  and  (S2,  E2, I2)  for  (2.1). 
Let  vs  =  Si  —  S2,  vE  =  Ei  —  E2  and  v1  =  Ii  —  h-  We  have  v5(x,0)  =  0  a.e.  in  [0,xmax], 
ve(x,9,  0)  =  0  a.e.  in  [0,xmax]  x  [0,oo)  and  vI(x,9,0)  =  0  a.e.  in  [0,  xmax]  x  [0,oo).  We  choose 
functions  <fs  £  Cl(Ds),  cj)E  £  Clr(DEI)  and  (j)1  £  C^r(DEI).  By  (2.2)  and  vs(x,  0)  =  0  a.e.  in 
[0,  xmax] ,  we  find  that  for  every  t  £  (0,  T) 

/“^max 

/  vs(x,t)4>s(x,t)dx 


rt 

=  /  <j)s(0 ,s)F(vs(-,s),vE(-,-,s),vI(-,-,s))ds 
Jo  . 


[4>s  (x,  s)  +  gS(x,  s)(fx(x,  s)  —  ms(x,  s)(fs(x,  s)]us'(x,  s)dxds 


'o,  Jo 

f*t  /“^ma  x 


(3.24) 


1  o±  Jo 

f*  t  (*Xm&x. 


vs(x,  s)<f>s(x,  s) 
S2(x,s)<fs(x,s) 


1 0  Jo 

/•^max  r°° 


'0  Jo 


\(y,  9)Ii(y,  9,  s)d9dy  dxds 
X(y,9)vI (y,9,  s)d0dy  dxds. 
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By  (2.3)  and  vE(x,9,  0)  =  0  a.e.  in  [0,  xmax]  x  [0,  oo),  we  have  that  for  every  t  G  (0 ,  T) 


max  /‘OO 


0  J  0, 


v^(x,  0 ,  t)cj)  (x,  9 ,  t)d0dx 

t  /‘lEmax  /*00 

[</>f  (x,  0,  s)  +  gE(x,  s)4>e ( x ,  0,  s)  +  4>E(x,  9,  s)]uE(x,  0,  s)d0dxds 


1 0  JO  JO 

f*t  /*^max  /*00 


[mE(x,  s)  +  /9E(0)]uE(x,  0,  s)(f>E(x,  9,  s)d6dxds 


io  Jo  Jo 

f*t  /^max 


+ 

+ 


'()  Jo 

f*t  /*Xmax 


vs(x,  s)cj)E(x,  0,  s ) 
-52 (x,  s)4>E(x ,  0,  s) 


/o  J 0 


^max  /‘OO 


UO  ^0 

3-max  /*00 


\(y,O)Ii(y,9,s)d0dy 

Hy,6)vI(y,0,s)d6dy 


dxds 

dxds. 


LJo  Jo 

By  (2.4)  and  id  (x,  0,0)  =  0  a.e.  in  [0,xmax]  x  [0,oo),  we  obtain  that  for  every  t  G  (0 ,T) 


3?max  /*  OO 


0  JO. 


v1  (x,0,t)(f)  (xiO^tydOdx 

t  ^max  /*00 

c/>7(x,  0,  s)  /  pE(9)vE(x,9,s)d9dxds 


io  Jp 

/*3?max  /*00 


/o 


+ 


'0  JO  J  0 

f*t  /*2)max  /*00 


[c/>f  (x,  9,  s )  +  </(x,  s)(f>x(x,  9 ,  s)  +  4>g(x,  9 ,  s)]x7(x,  0,  s)d9dxds 


m1  (9)vI (x,  0,  s)</d  (x,  0,  s)d9dxds. 


0  JO  Jo 


We  choose  a  function  <(r  G  C'1(D5)  that  satisfies  the  problem 

</>f  (x,  s)  +  gS(x,  s)4>E(x,  s)  =  0,  0  <  s  <  t,  0  <  x  <  xE 

41  (®max>  s)  —  0,  0  <C  S  <C  t, 

(t>S{x,t)  =  Ts(x),  0  <  x  <  xmax, 


(3.25) 


(3.26) 


(3.27) 


where  Y5  G  Cq°((0,  xmax))  and  —  1  <  Ts  <  1.  We  also  choose  a  function  <f>E  G  Cqt{Dei )  that 
satisfies  the  problem 


<l>f  (x,  0,  s)  +  gE{x,  s)(j)E (x,  0,  s)  +  (j)E (x,  0,  s)  =  0,  (x,  0,  s)  G  (0,  xmax)  x  (0,n)  x  (0 ,t), 

4>E{x max,  0,  s)  =  0,  0  <  S  <  t,  0  <  0  <  n, 

4>E(x,  n,s)  =  0,  0  <  S  <  t,  0  <  X  <  Xmax, 

4>E(x,  9,  t)  =  Te(x,  0),  0  <  x  <  xmax,  0  <  0  <  n, 


(3.28) 


where  Ye  G  Co°((0,  xmax)  x  (0,  n))  and  —1  <  TE  <  1.  We  further  choose  a  function  (j)1  G  Cqt(Dei ) 
that  satisfies  the  problem 


<j>l(x,  0,  s)  +  /(x,  s)^(x,  0,  s)  +  (f)g(x,  9,  s)  =  0,  (x,  0,  s)  G  (0,  xmax)  x  (0,  n)  x  (0,  t), 

0^(j'max;  9,  s)  =  0,  0  <  s  <  t,  0  <  0  <  n,  (3  29) 

07(X,  U,  S)  =  0,  0  <  S  <  t,  0  <  X  <  Xmax, 

07(X,  0,  t)  =  YJ(X,  0),  0  <  X  <  Xmax,  0  <  0  <  n, 

where  Yy  G  Co°((0,  xmax)  x  (0 ,n))  and  —  1  <  TJ  <  1.  The  existence  of  such  cf)S(x,s),  4>E(x,9,s), 
and  cj)I( x ,  0,  s)  is  guaranteed  since  all  the  above  problems  are  linear  with  local  boundary  conditions. 
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The  boundary  and  initial  conditions  of  the  above  three  problems  imply  that  —1  <  (f>s ,  (pE ,  J)1  <  1. 
Substituting  4>s  that  satisfies  (3.27)  into  (3.24),  we  find 


<  iAs 


vs  (x,  t)Th(x)dx 

f*t  /*3?max 


|u's(a:,  s)\dxds  +  vS'E 


t  COO 


0  Jo 


+U 


S,I 


t  r^max  coo 


\vE(x,  0,  s)\d0dxds  (3.30) 


o  Jo  Jo 


\vJ (x,  0 ,  s)\d0dxds, 


o  Jo  Jo 


where  uS’S  =  ||/3'S,||0o  + l|w5|loo  + 


oo  W  IK' I  loo 


sup 


f*^max  pOO 


h (y,9,s)d6dy  L  vS’E  =  ||/?B||oo  and  uSJ  = 


|oo  ~r  IK'Hoo  SUp 

se[0,T]  iJo 


se[o,T]  l  Jo  Jo 

S2{y,  s)dy\.  Substituting  cf>E  that  satisfies  (3.28)  into  (3.25),  we  have 


>0  Jo 

<VE’S 


vE(x,  0,  t) TE(x,  0)d0dx 


t  rx-n 


t  /“^max  POO 


0  JO 


pt  pX max  pOO 

+uE’1  /  /  /  \vI(x,0,s)\dOdxds, 

Jo  Jo  Jo 


vs(x,  s)\dxds  +  uE,E  II  I  \vE(x,0,s)\d0dxds 

Jo  Jo  Jo 


(3.31) 


where  the  constants  v  E-S  =  Mlco  sup 


r^-max  pOO 


Ii(y,0,s)d0dy'>,  vE’E  =  ||mii||00  +  \\pE  ||oo  and 


vE’T  = 


oo  sup 
se[o,T]  l Jo 


se[o,T]  Uo  Jo 

S2(y,  s)dy\ .  Substituting  < j)1  that  satisfies  (3.29)  into  (3.26),  we  obtain 


f^max  pH 


v\x ,  0,  t) T1  (x,  0)d0dx 


t  pX max  POO 


1 0  Jo 

<uI,E  If  I  \vE(x,  0,  s)\d0dxds  +  v1'1 

Jo  Jo  Jo 


t  PX  max  POO 


\vk(x,  0 ,  s)\d0dxds , 


(3.32) 


0  JO  JO 


where  the  constants  uI,E  =  ||pE||oo  and  u1'1  =  ||mJ||oo. 

Since  (3.30)  holds  for  every  Y5,  we  can  for  each  fixed  t  choose  a  sequence  {Y5,fe}^=1  converging 


to  TS ,  where 


(  1  if  vs(x,  t)  >  0 

TS(x)=i  0  ifus(x,f)=0  . 

[  —1  if  vs(x,  t)  <  0 

Substituting  TS'k  into  (3.30)  and  using  the  dominated  convergence  theorem,  we  find 


|u‘s'(a;,  t)\dx 


t  rx  n 


<vS'S  /  /  \vs(x,s)\dxds  +  us,E 


t  PX  max  COO 


0  Jo 


+v 


S,I 


t  C3-max  COO 


\vE(x,  0,  s)\d6dxds  (3.33) 


o  Jo  Jo 


\vJ (x,  0,  s)\d0dxds. 


o  Jo  Jo 
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We  also  note  that  (3.31)  holds  for  every  TE ,  hence  we  can  for  each  fixed  t  choose  a  sequence 
{TE,k}(^L1  converging  to  TE ,  where 


rE(x,e) 


1  if  vE(x ,  9,t)  >  0 

0  if  vE(x,  6,t)  =  0 

—  1  if  vE(x,  6,t)  <  0 


Substituting  such  TE'k  into  (3.31),  using  the  dominated  convergence  theorem,  and  using  the  fact 
that  the  constants  uE’S ,  isE'E  and  vE’1  are  independent  of  n,  we  find 


nCXD 

<vE’S 
.EJ 


\vE(x,  0,  t)\dOdx 


|us(a:,  s)\dxds  +  uE'E 


t  /“^max  P  °C 


0  JO 


t  /“^max  pOO 


\vE(x,  9,  s)\dOdxds  (3.34) 


o  Jo  Jo 


+v 

Similarly,  we  can  find 


| v1  (x,  6,  s)\d9dxds. 


o  Jo  Jo 


f*X max  POO 


'0  JO 


<  V 


IE 


\v\x,  6,  t)\dOdx 

(*t  pX max  POO 


t  PXrnax.  POO 


\vE(x,  9,  s)\d9dxds  +  v1’1 


o  Jo  Jo 


|?/(x,  6,  s)\d8dxds. 


o  Jo  Jo 


/•OO 

By  the  fact  /  e~edQ  =  1  and  (3.33)-(3.35),  we  have 

Jo 

[e~9\vs(x,  f)|  +  \vE(x,  9,t) |  +  [v1  (x,  9 ,  t)\\d9dx 

) 

[e~9\vs(x,  s)|  +  \vE(x,  9,  s)|  +  \v\x ,  9,  s)\]d9dxds, 


(*X max  POO 


/  0  Jo 


t  PXrna,x  POO 


(3.35) 


<  v 

1 0  Jo  Jo 

where  v  =  max{i/s,s  +  vE,S ,  vS,E  +  uE,E  +  uI,E,  vS)I  +  vE'1  +  v1’1}-  By  Gronwall’s  inequality,  we 
thus  have 


f^max  pOO 


[e  0|i>s(x,i)|  +  \vE(x,6,t) |  +  | v1  (x,9,t)\\d0dx  =  0. 


Jo  Jo 

n  Q  TP  T  TP  T 

Hence,  we  have  v  =  0  a.e  in  D  ,  v  =  0  and  v  =  0  a.e.  in  D  . 


□ 


4  Monotone  Approximation  and  Existence 

Suppose  that  (S°,E°,I°)  and  (S°,E°,I°)  are  a  pair  of  nonnegative  upper  and  lower  solution  of 
(2.1),  respectively,  then  by  Theorem  3.2  we  have 

S°(x,t )  <  S°(x,t )  a.e.  in  Ds , 

E°(x,9,t )  <  E°(x,9,t )  a.e.  in  .DE/,  (4.1) 

I°(x,9,t)  <  I°(x,9,t)  a.e.  in 
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We  set  up  two  sequences  {<Sfc,  Ek,  I_k}kLi  and  {S"',  E*,  by  the  following  procedure: 


-pzk  -= fc 


'S’i  (M)  +  {gs(x,t)S  (x,t))x  =  -S  (x,t) 


\ 


■7=7/2 — 1  /  A  ~r=;k— 1. 


ms(x,  t )  + 


-p/c— 1  . 


s-vk 


S  (x,  0)  =  Su(x), 


A  (y,e)Ik-1(y,9,t)d9dy 


0  ■'  0 


(4.2) 


S%( x,  t )  +  (gS(x,  i)<Sfc(x,  t))a;  =  —S_k(x,  t)  ms (x,  t )  + 
/(0,t)5fc(0,t)  =  F(5fc-1(At),^-1(v,t),/fc-1(-,-,t)) 
Sfc(x,0)  =  S°(x), 


Hy,0)Ik  1(y,0,t)d6dy 


0  Jo 


(4.3) 


Et  (x,9,t)  +  (gE(x,  t)E  (x,9,t))x  +  Ee(x,9,t)  =  ~[mE(x,t)  +  pE{9)}E  (. x,9,t ) 

_ u  —k-l  fXmaxroo  _ 

£  (x,0,t)  =  S'  (x,i)  /  /  X(y,9)I  (y,9,t)d9dy 

_  _  40  4o 

E  (0,9,t)  =  0,  E(s,0,O)  =  EP{x,0), 


(4.4) 


Ek(x,  9,  t)  +  (gE(x,t)Ek(x,9,t))x  +  E%(x,9,t)  =  ~[mE(x,t)  +  pE (9)]Ek(x,9,t) 


E*(x,0,t)  =  Sk~\x,t) 


A  (y,9)Ik-\y,9,t)d9dy 


,  -  0  J  0 

J7fc(0,  9,  t)  =  0,  Ek(x,  9 ,  0)  =  E°(x,  9) 


(4.5) 


It(x,9,t )  +  (gI(x,t)Ik(x,9,t))x  +  Ig(x,9,t)  =  - mI(9)Ik(x,9,t ) 

7fc(x,0,f)=  pE(9)Ek-\x,9,t)d9  (4.6) 

40  _ 

r  (0,  0,  t)  =  0,  r(x,0,O)  =  7°(x,0), 

It  (»,  0,  i)  +  (9J(ar,  t)Ifc(x,  9,  t))x  +  Ig(x,  9 ,  t)  =  -m7(60/fc(x,  9,  t) 

/•OO 

lk(x,0,t)  =  pE(9)Ek~1(x,  9,  t)d9  (4-7) 

/fc(O,0,t)  =  O,  Ik{x,9,0)  =  7°(x,0). 

The  existence  of  the  solutions  to  the  above  problems  (4.2)— (4.7)  follows  from  standard  results,  given 
the  fact  that  they  are  all  linear  problems  with  local  boundary  conditions  for  Sinko-Streifer  type 
systems. 

We  first  show  that 

S°(x,t )  <  S}(x,t)  a.e.  in  Ds , 

E°(x,9,t)  <  ^(x,  9,  t)  a.e.  in  DEI ,  (4.8) 

I°(x,9,t )  <l}(x,  9,t)  a.e.  in  DEI . 
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Let  us  =  S°  —  S1,  uE  =  E°  —  E1  and  u1  =  1°  —  I1,  then  by  (4.3),  and  using  the  fact  that  (S°,  E°,  1°) 
is  a  lower  solution  of  (2.1),  we  have 


f^max 


< 


us(x,  t)ijjs(x,  t)dx 

t  (*%  max 

[ifjg  (x,  s )  +  gs(x,  s)ipE(x,  s)  —  ms(x,  s)'tps(x,  s)]u5'(a:,  s)dxds 

t  r^max 


0  J, 0 


TTmax  POO 


Ub(x,s)llj  (x,s) 


0  JO 


A (y,  9)1° (y,  9,  s)d9  dy  dx  ds. 


o  Jo 


Hence,  uS(x,t)  satisfies  (3.22)  with  As(x,t)  =  —ms(x,t)  ~ 


3? max  /*00 


A (y,9)I°(y,9,t)d9dy.  Thus, 


o  Jo 


S°(x,t)  <  S}(x,  t).  By  (4.5),  and  using  the  fact  that  (S° ,  E° ,  7°)  is  a  lower  solution  of  (2.1),  we 
have 


31  max  /*00 


uE(x,  9,  t)^E(x,  9,  t)d9dx 


0  JO, 


t  /*31ma x  /*00 


< 


[ipf(x,  9,  s )  +  gE(x ,  s)i/je(x,  9,  s )  +  ipo(x,  s)]uE(x,  9,  s)d9dxds 


/ o  Jg  Jo 

c*  rx max  roc 


[mr  (x,  s)  +  p  ( 9)\yr  (x,  0,  s)%j)  (x,  0,  s)dOdxds. 

10  Jo  Jo 

Hence,  uE(x,9,t )  satisfies  (3.23)  with  AE(x,9,t)  =  — mE(x,t )  —  pE (9).  Thus,  E°(x,9,t )  < 
H1( x,9,t).  By  (4.7),  and  using  the  fact  that  (5° ,  E° ,  7°)  is  a  lower  solution  of  (2.1),  we  have 


31  max  pOC 


to  Jo, 

< 


u7(x,  9,  t)^(x,  9,  t)d9dx 

t  /*  3J  max  /*00 

[V’f  (a;,  0,  s)  +  g1  (x,  s)t/j^(x,  9,  s )  +  ipg(x,  9,  s)]uJ (x,  9,  s)d9dxds 


10  Jo  Jo 

rt  f*X  max  /*00 


mI(9)uI(x,  9,  s )i/j1  (x,  9 ,  s)d9dxds. 


o  Jo  Jo 


Hence,  I°(x,9,t )  satisfies  (3.23)  with  AE(x,  9,  t)  =  —mI(9).  Thus,  I°(x,9,t )  <  l}(x,9A)-  Similarly, 
we  can  show  that 


"7x0 


"jxl 


5  (x,f)  >  S'  (x,  t)  a.e.  in  IA  , 


x=0 


Bu(x,0,f)  >  El(x,0,t)  a.e.  in  D 


EI 


(4.9) 


v0 


/u(x,#,t)  >  E(x,9,t)  a.e.  in  ZA^1. 


-fL 


By  (4.2)-(4.9),  we  can  see  that  (S^,  E^ ,  E)  and  {S} ,  E1 , 11)  are  a  pair  of  upper  and  lower  solutions 
of  (2.1),  respectively.  Hence,  by  Theorem  3.2  we  have 


S}(x,t)  <  S1(x,t)  a.e.  in  Ds , 

Ei(x,9,t)  <  E1(x,9,t)  a.e.  in  IA57,  (4.10) 

l}(x,Q,t)  <  /1(x,  0,t)  a.e.  in  ZA^1. 
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We  then  assume  that  for  some  k,  ( Sk,Ek,Ik )  and  ( Sk,Ek,Ik )  are  a  pair  of  upper  and  lower 
solutions  of  (2.1),  respectively.  Proceeding  analogously,  we  can  first  show  that 

S_k(x,t )  <  *Sfc+1(x,i)  and  Sk(x,t)  >  Sk+1(x,t )  a.e.  in  Ds , 

Ek(x,8,t )  <  Ek+1(x,  8,  t)  and  Ek(x,8,t)  >  Ek+1(x,8,t )  a.e.  in  DEI ,  (4-11) 

I_k(x,6,t )  <  Ik+1(x,8,t)  and  Ik(x,8,t )  >  Ik+1(x,8,t)  a.e.  in  DEI , 

and  then  by  (4.2)-(4.7)  and  (4.11),  we  can  claim  that  (Sk+1,  Ek+1,  Ik+1)  and  (Sk+1,  Ek+1,  fk+1)  are 
a  pair  of  upper  and  lower  solution  of  (2.1),  respectively.  Thus,  we  obtain  two  monotone  sequences 
{Sk,Ek,Ik}  and  {Sk,EkJk}  which  satisfy 

S°  <  S1  <  ■  ■  ■  <  Sk  <  ■  ■  ■  <  Sk  <  ■  ■  ■  <  S1  <  S°  a.e.  in  Ds , 

E°  <  E1  <  ■  ■  ■  <  Ek  <  ■  ■  ■  <  Ek  <  ■  ■  ■  <  El  <  E°  a.e.  in  DEI , 

Z°  <  I1  <  •  •  •  <  Ik  <  ■  ■  ■  <  Ik  <  ■  ■  ■  <  I1  <  7°  a.e.  in  DEI . 


From  the  nronotonicity  of  the  sequences  { , EK , I_; }  and  {S  ,E,I  }  it  follows  that  there  exist 
functions  (S,g,I)  and  (S,E,I)  & uch  that  (gk,EkJk)  ->  (S,E,I)  and  (Sk,Ek,Ik)  -»•  (S,E,1). 
Clearly,  S_  <  S  a.e.  in  Ds ,  E<  E  and  /  <  I  a.e.  in  DE1 .  On  the  other  hand,  by  the  dominated 
convergence  theorem,  we  know  that  (S_,E,E)  and  ( S,E,I )  are  also  an  upper  and  lower  solutions 
of  (2.1),  respectively.  Hence,  by  Theorem  3.2,  we  have  5  >  S  a.e.  in  Ds ,  E>E  and  I  >  I  a.e. 
in  DE1  .  Thus,  S  =  S  a.e.  in  Ds ,  E  =  E  and  /  =  I  a.e.  in  DE1 .  Defining  this  common  limit 
by  ( S,E,I ),  then  we  see  that  it  satisfies  (2.2)— (2.4)  for  every  nonnegative  functions  G  C|1(D5'), 
i\)E  G  C'lr{DEI)  and  ^  G  Cq  r(DEI).  By  the  property  of  modifiers  and  the  dominated  convergence 
theorem,  we  can  easily  show  that  the  limit  satisfies  (2.2)— (2.4)  for  every  function  (j)S  G  Cll(D5'), 
4>e  G  Cqt(Dei)  and  (b1  G  Cqt{Dei).  Hence,  ( S,E,I )  is  the  solution  of  (2.1). 

Based  on  the  above  discussion  and  Theorem  3.4,  we  have  the  following  existence-uniqueness 
result. 


Theorem  4.1.  Suppose  that  all  the  assumptions  (H1)-(H5)  hold.  We  assume  that  (S°,E°,I°) 
and  (S°,E°,I0)  are  a  pair  of  nonnegative  upper  and  lower  solution  of  (2.1),  respectively.  Then, 
there  exist  monotone  sequences  {5fc ,  Ek ,  fk}  and  {5  ,E,I}  which  converge  to  the  unique  solution 
(S,E,I)  of  (2.1). 

Remark  4.2.  For  the  initial  data  S°(x)  =  r]S  exp(— asx),  E°(x,8)  =  r]E  exp(— aEx)  exp(— 8)  and 
I°(x,  8)  =  rj1  exp(— oJx)  exp(— 9)  with  some  positive  constants  as ,  aE ,  a1,  r/s ,  r/E  and  g1 ,  we  can 
construct  the  nonnegative  lower  and  upper  solution  as  follows: 


(S(x,t),  E(x,8,t),  I(x,e,t))  =  (0,0,0), 
S(x,  t )  =  as  exp (bst.)  exp(— dsx), 

E(x,8,t)  =  - 7^aE  ex.p(bEt)  exp(— dEx), 

1  +  8- 

I(x,8,t )  =  - -tg/  exp(6Jf)  exp(— dkx), 

1  +  8Z 
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q  zp  T  ^  E1  T  Q  fT  T 

with  positive  constants  a  ,  a  ,  a  ,  6,6,6,  d,d  and  d  which  are  to  be  determined.  First  we 
choose  positive  constants  bE  and  b1  large  enough  such  that 

6s  >  1  +  max  I  of' (a;,  t)|  +  a7'  max  gE(x,  t), 

(: x,t)eD  (x,t)£D 

b1  >  1  +  max  |g7(x,  t)|  +  a7  max  g\x,  t), 

(x,t)£D  (x,t)£D 

where  D  =  [0, 1]  x  [0,T].  Fixing  these  constants  bE  and  b1 ,  then  we  choose  6s  sufficiently  large  so 
that 

bs  >  max  <  max  \g~(x,  t)|  +  as  max  gs(x,  t ),  6®,  6 1 
[(xfytD  (x,t)£D 

The  positive  constant  aE  is  chosen  so  that  aE  >  r]E ,  and  a 1  is  chosen  so  that 

a1  >  max|r//c^||pE||0Oexp(6ET)aE|  . 

The  positive  constant  as  is  then  chosen  to  satisfy  as  >  {gs ,  aE ,  a1}.  We  then  fix  imai  sufficiently 
small  so  that 


<  min  <  1 


2  minte[0ir]  gs(0,  t )  4  minte[0)T]  gs(0,  t ) 


L  9||/3s||oo 

4mintg[0iy]  gs(0,  t) 
97r||/37||oo 


97t||/3^ 


4cr 


’  37r|| AHooa^a7 


exp (—6s  -  6/)T) 


where  the  L ^  norm  of  (5s  is  defined  on  (0, 1),  the  L ^  norms  of  f3E ,  /31  and  A  are  all  defined  on 

(0,1)  x  (0,oo).  Note  that  since  hm  - —  — mdX ~  =  xmax,  we  can  choose  ds  sufficiently 

small  so  that 

ds  s  oS  md  1  -  exP(-dsxmax)  s  3X_ 

Similarly,  dE  is  chosen  sufficiently  small  enough  so  that 

and  dE  <  min{a:E,  ds}. 


1  -  exp (-dExmax)  ^  3 

dE  ~  2Xmax 


Fixing  these  constants  for  ds  and  dE ,  then  we  choose  d1  such  that 


1  -  exp(-d7xmax) 
d1 


<  -xmax  and  cr  <  min{or ,  dE}. 


With  all  these  positive  constants  fixed  as  above,  we  can  easily  show  that  (S_(x,  t),E(x,  9,  t),I(x,  6 ,  t)) 
and  (S(x,t),  E(x,9,t),I(x,6,t))  are  a  pair  of  nonnegative  lower  and  upper  solutions  of  (2.1),  re¬ 
spectively. 
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5  Concluding  Remarks 


In  the  above  presentation  we  developed  a  comparison  principle  and  constructed  monotone  sequences 
to  establish  the  existence  and  uniqueness  of  solutions  for  a  class  of  SEI  epidemic  models  in  which 
individuals  in  the  exposed  state  and  infectious  state  are  structured  by  size  and  class  age  (or  residency 
time).  A  special  case  of  the  general  class  is  the  biomass/ viral  infection  model  of  [6].  We  believe 
that  this  method  can  be  applied  to  a  more  generalized  version  of  this  model  in  which  mortality 
rates  and  birth  rates  depend  on  the  total  population  or  total  biomass  of  each  state;  this  would  yield 
results  under  even  more  reasonable  and  realistic  assumptions.  Note  that  we  defined  coupled  upper 
and  lower  solutions  for  the  equation  in  the  susceptible  state  because  of  the  nonlinear  term  in  this 
equation.  If  we  wish  to  apply  this  method  to  a  more  general  case,  we  may  need  to  define  totally 
coupled  upper  and  lower  solutions  for  the  equations  of  all  the  states.  Because  of  the  complexity  of 
our  model,  there  are  other  issues  we  have  not  considered  in  this  paper  such  as  calculating  the  basic 
reproduction  ratio  to  determine  if  the  disease  persists  or  dies  out.  Nor  have  we  considered  control 
of  the  population  environment  to  yield  maximum  (for  the  case  of  the  biomass/ vaccine  production 
problem)  or  minimum  (for  the  case  of  epidemic  prevention  among  healthy  populations)  infected 
biomass.  These  and  other  issues  are  the  subject  of  current  efforts. 

There  are  other  methods  available  in  the  literature  that  can  be  used  to  prove  existence  and 
uniqueness  of  solutions  to  problems  such  as  (2.1).  We  have  discussed  some  of  these  in  the  intro¬ 
ductory  section.  Weak  solutions  and  semigroup  theory  have  been  used  in  [7,  8]  to  prove  existence 
and  uniqueness  for  a  classical  size  structured  population  model  by  computing  the  existence  of  the 
solution  semigroup  in  an  extrapolation  space.  In  [9]  a  classical  size  structured  population  model 
that  includes  non-observable  characteristics  responsible  for  variations  in  growth  rates  for  individ¬ 
uals  of  the  same  size  is  investigated.  Well-posedness  of  this  model  and  a  nonlinear  perturbation 
of  it  is  proved.  Our  future  endeavors  will  include  investigation  of  the  use  of  weak  formulations 
and  nonlinear  semigroup  theory  to  establish  existence,  uniqueness  and  continuous  dependence  of 
solutions  to  problems  such  as  (2.1).  A  sensitivity  analysis  framework  is  also  under  development. 
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